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Abstract 


NASA is developing a systems biology approach to improve the assessment of health 
risks associated with space radiation. The primary toxic and mutagenic lesion 
following radiation exposure is the DNA double strand break (DSB), thus a model 
incorporating proteins and pathways important in response and repair of this lesion is 
critical. One key protein heterodimer for systems models of radiation effects is the 
Ku70/80 complex. The Ku70/80 complex is important in the initial binding of DSB 
ends following DNA damage, and is a component of nonhomologous end joining 
repair, the primary pathway for DSB repair in mammalian cells. The SAP domain of 
Ku70 (residues 556-609), contains an a helix-extended strand-helix motif and similar 
motifs have been found in other nucleic acid-binding proteins critical for DNA repair. 
However, the exact mechanism of damage recognition and substrate specificity for 
the Ku heterodimer remains unclear in part due to the absence of a high-resolution 
structure of the SAP/DNA complex. 

We performed a series of molecular dynamics (MD) simulations on a system with the 
SAP domain of Ku70 and a 10 base pairs DNA duplex. Large-scale conformational 
changes were observed and some putative binding modes were suggested based on 
energetic analysis. These modes are consistent with previous experimental 
investigations. In addition, the results indicate that cooperation of SAP with other 
domains of Ku70/80 is necessary to explain the high affinity of binding as observed in 
experiments. 



Introduction 
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About 5x1 0 5 Ku molecules per nucleus. 

Ku70 is mainly responsible for DNA 
binding. 

Two DNA-binding domains: N terminus 
before residue 440, C terminus (SAP 
domain) 



DSB repair by nonhomologous 
end joining (NHEJ) 


Helix-extend strand-helix 
structure of SAP domain 


Material and Methods 


Resources: 

• Computational programs package: AMBER 9.0 

• Workhorse: Radiation Biology Laboratory Beowulf Cluster (16X4 
AMD Opteron x86-64) 

• SAP domain: NMR structure (PDB accession code 1 JJR) 

• DNA 10-mer: d(GCGTTAACGC) 2 , generated by NUCGEN of 
AMBER 

MD protocol: 

1. Putative binding modes search: implicit solvent (GB2), multiple 
alternative constraint/restraint simulation annealing (SA) (250ps) 
and semi-flexible MD (250ps) 

2. Binding affinity calculation: 5ns unconstraint explicit solvent (TIP3P 
water), counterions (Na + ), 300K, 1 atm, SHAKE, PME, MM_PBSA 
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Force Fields and Solvent 

Models 


Performance of three force fields on SAP 
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Performance of three implicit solvent on SAP 



RMSD: root mean square deviation (A) of heavy atoms of backbone. 
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How about DNA? 


Unconstraint simulation of SAP with ff99SB and TIP3P water Unconstraint simulation of DNA with GB2 implicit solvent 




DNA is more flexible than SAP 



Starting Structures (manually 

docked) 


kuc3 


kuc4 


kuc5 


Chemical shift changes upon 
complex formation ( in yellow) and 
conserved Arg or Lys residues. 
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Potential energy (kcaJ/mol) 


Constraint/Restraint SA 


Potential energy variation along the path of SA 
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Annealed Structures (most visited 

conformations) 



> Loops rather than helices directly interact with the 
minor groove of DNA. 

>DNA experiences apparent bending upon binding. 
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Explicit Solvent Simulation 


RMSD (A) of heavy atoms of the complex backbone Separate RMSD (A) of SAP and DNA 
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>Complex is stabilized after 2 ns equilibration. 

>There is apparent adjustment when the complex is immersed into 
explicit solvent, but the binding mode is well maintained. 


Binding Affinity Calculation 
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Binding Free Energy Components 
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a Contribution (SAP-DNA)-contribution (SAP)-contribution (DNA). 
b kcal/mol. Average over 150 (15 in the case of entropy contributions) snapshots. 
c Standard error of mean values. 


Conclusions and Future Work 


SAP domain of Ku70 is unlikely to be involved in DSBs end 
recognition, but shows apparent affinity to DNA minor groove. 

The loops rather than the helices of the SAP domain are more 
active in binding DNA, and in stable complex, the DNA substrate 
is bent while the structure of SAP is well preserved. 

The converged stable complex found by implicit solvent method 
are verified by explicit solvent computation, indicating a 
reasonably effective way to search the possible binding modes 
between macromolecules. 

Different sequences of DNA duplex will be investigated to 
understand the mechanism of pause sites when Ku translocating 
along DNA. 



